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Abstract 

We present a scheme to realize the chiral topological excitonic insulator in semiconductor het- 
erostructures which can be experimentally fabricated with a coupled quantum well adjacent to two 
ferromagnetic insulating films. The different mean-field chiral topological orders, which are due to 
the change in the directions of the magnetization of the ferromagnetic films, can be characterized 
by the TKNN numbers in the bulk system as well as by the winding numbers of the gapless states 
in the edged system. Furthermore, we propose an experimental scheme to detect the emergence 
of the chiral gapless edge state and distinguish different chiral topological orders by measuring the 
thermal conductance. 
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I. INTRODUCTION 



The search for new phases of quantum matter is one of the essential topics in condensed- 
matter physics. Chiral topological band insulators (TBIs) are such a type that has been 
attracting a lot of interest both theoretically and experimentally Although like trivial 
insulators in the sense that TBIs have a band gap in the bulk, they are fundamentally 
distinguished from trivial ones by their having gapless modes on the boundaries. These 
gapless modes are robust under perturbations and cannot be gapped without going through 
a quantum phase transition. In the case of time reversal symmetry (TRS) breaking, a well- 
known TBI system is the Haldane's model which is a minimal model to illustrate quantum 
anomalous Hall effect (QAHE) [1]. The QAHE topological phase is characterized by the 
(Thouless, Kohmoto, Nightingale, and Nijs) TKNN number [2] of the first Chern class of 
a U(l) principal fiber bundle on a torus in the bulk system or the winding number of 

n n 

Halperin's edge-state theory |3|, |4| on the boundary of the system. The coherence of the 
two different kinds of numbers is guaranteed by the bulk-edge correspondence. Since the 
rigorously prerequisite magnetic field in Haldane's model is difficult to realize in experiment, 
recently, there are some new proposals [5] to realize QAHE based on single-particle picture. 
In analogy with QAHE in single-particle picture, the superconductors in TRS-broken 



e bulk gap opened by 
6j. The edge states of 



{Px + iPy) weak pairing state in two dimensions with a fully stab 
electron-electron interaction can also have chiral topological order 
the chiral superconductor have half of the degrees of freedom compared to QAHE states due 
to the particle-hole symmetry (PHS) and are called Majorana edge sates. In the spirit of 
analogy with superconductor, a natural and important issue is how to get chiral topological 
excitonic insulator (TEI), which is addressed in this paper. 

In this paper, we consider an independently gated double-quantum-well structure sepa- 
rated by a spacer as shown in Fig. 1(a). The ferromagnetic insulating films are introduced to 
break TRS by inducing an effective Zeeman splitting in the two dimensional electron (hole) 
gas [2DE(H)G]. The magnetization is perpendicular to the two-dimensional layer. Note that 
the orbital effect of the ferromagnetic films to the 2DE(H)G can be neglected due to the 
local exchange interaction on the interface. An electron-hole fluid is created by modulating 
the voltages so that the Zeeman-split upper branch of the heavy-hole bands in 2DHG layer 
can move above the Zeeman-split lower branch of the electron bands in 2DEG layer. This 
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FIG. 1: (Color online) (a) Schematic structure of semiconductor quantum- wells system that holds 
chiral TEL The external gates (V ge (h)) can independently tune the chemical potential [i e {K\ to 
obtain the electron and heavy-hole layer. The ferromagnetic insulating films support effective 
exchange fields (V e ,Vh). (b) The energy spectrum of the electron/hole bilayer system near the 
Fermi energy Ep ~ 0. Here the solid lines denote non-interacting single-particle energy spectrum 
£> e (7i)±> while the dashed lines denote the exciton energy spectrum with an obvious mean-field gap 
opened. We take t e = th = 1, /x e = = ~ 4, (V e , Vh) = (1, 1), a = 0.5, and (3 = 0.5. 

procedure results in spatially separated but strongly interacting electron and hole fluids if 
the two layers are close enough. The external electric field produced by the bias voltages 
and the intrinsic electric field due to doping in the process of fabricating quantum wells can 
enhance the structural inversion asymmetry and induce the tunable Rashba spin-orbit (SO) 



interaction 



7 



10(. Recall that due to the strong SO coupling and non-centrosymmeytic 
property |rjj, the broken parity of the order parameter is the prerequisite condition of the 
chiral superconductor. In analogy with chiral superconductor, we demonstrate that the chi- 
ral TEI can occur when the Rashba SO interaction is strong enough with respect to the 
amplitude of the excitonic order parameter (EOP). Generally the Rashba SO interaction 
strength can be influenced by carrier density, gated voltages, material of quantum wells, etc 



12Hl4| . In InAs heterostructures, for instance, this quantity can be electrically tuned to be 



as large as a~50 meV A 15]. Due to the missing PHS, the edge states of chiral TEI in the 
present system are not Majorana fermions. The implication of these gapless edge states for 
experimental observations is also discussed in this paper. 



3 



II. MODEL AND HAMILTON! AN 



We start with an effective electron/hole semiconductor bilayer system confined in the 
x-y plane. Here for the hole layer only the heavy-hole bands are occupied as in typical 
experiments, while the light-hole bands are empty and are therefore not taken into account 
in our model. The resultant tight-binding Hamiltonian for the Rashba spin-orbit coupled 
semiconductor bilayer system is %=X!p(^^i+^R ) )+^!n7 h) =^o±% i t7 h): 

<i,j>,(T j,a,a' 
3 

~ E( e lt e i+ 5 ^ ~ e h e 3-8y,i)\ + Hx -> 

3 

3 

3 

+ 3(1 - i] S ^{h] 1[ h j -8 X+ 8y A - h^hj+sx-sy,^) (1) 

3 

+ 3(1 + i)y~]{h\,ithj-5x-6y,i). - h] :t h j+Sx+ sy^) 
j 

3 

j+sy,n - h)^hj-sy^)\ +H.c, 

3 

n^ H) --- V U {eh) (d)e ] h j h- ,e- 

n int — 2 2-^ hi \ a ) (i ia"-ja''h<7' (i w 
i,3,a,a' 

Here t p denotes the nearest-neighbor hopping amplitude while ji p represents the chemical 
potential in electron (p—e) or heavy-hole {p—h) layer. s z is the z-component of the Pauli 
matrices and V p r p is the effective Zeeman splitting (r e =l for electron layer and 7^=— 1 for 
hole layer). pj a is the fermion annihilation operator at lattice site j with spin ±1/2 (t, i) 
for p=e and spin ±3/2 (fT,JJ-) for p—h. a (j3) is the Rashba SO interaction strength in the 
electron (heavy-hole) layer. 5x (8y) is the square lattice spacing along the x (y) direction. 

In the interaction term, U^ J h \d)—e 2 j 'e\j Ri >e — Rj,h +d 2 , where e is the dielectric constant 
of the spacer and d is the interlayer distance. We only consider the interaction correla- 



tive to exciton formation and ignore the electron-hole exchange interaction. The lattice 
Hamiltonian can be transformed into the momentum space with the Fourier transformation 

( e k,a'h,*) = 1 /^T,i eik ' Rl ( e i,v, h i,«)- The reSult reads 

k,a,a',p 

7~^Fi(k) — y~^a:(sin k x — % sin k y )el + H.c, 



nj{e-h) J_ V- ^(efc)/ x f fe t ^ g ^ 



(2) 



k,k' ,q,a,a' 

where JJ <yeh \q) = '^-e'~ qd , C?^=— 2t p (cos A^+cos A^), afc=2(3 cos A^ sin A; x — sin fc x cos A; x — 2 sin fc x ), 
and bk=2(— 3 cos fc x sin Aj^+sin k y cos /c y +2 sin In the above Hamiltonian, the interlayer 
tunneling is neglected, because the insulating spacer can supply a high barrier to stop the 
direct interlayer hopping. We also neglect the intralayer electron-electron and hole-hole 
interactions, since they are expected to renormalize the single-particle energy of each 
layer and have no essential influence on the topological properties of the system. In the 
mean-field approximation, the above Hamiltonian can be written as 



k,cr,cr' ,p 



+ £<-J£(^*(*>1^+h.c.) 



(3) 



kcrcr' 



+ 



-£ 



kqcrcr' 



Kcn{k)KAk-q) 
U( eh 1{q) 



where the EOPs are defined as 



(4) 



In the Nambu notation with combined e-h field operator basis ib=\£i+ er, h - h - 1 T , the 
mean- field Hamiltonian is expressed as T-Lmf—^ ' Hmf^^ const with 

/i e + \/ e s 2 A(fc) 



MF 



e<?> 



At(jfe) 



(5) 
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where Y l {p l=±C > {p) -.I±H { E ) (±k) and 

±k ±k n y ' 

AW = -\ 

with A aa ,(k) defined in Eq. (|4]). 

The complex EOPs A acr /(k) can be self-consistently obtained from exact numerical cal- 
culation of Eqs. ([3]) and (J4]) with respect to minimizing the ground state energy. In our 
numerical calculation of A aa i(k), we set the lattice size 81x81 and take t e =t h =l, fi e =fi h =—4:, 
a=0.5, and (3=0.5. There are four different kinds of choices for the perpendicular magnetiza- 
tion in the two magnetic films adjacent to the bilayer system. For the parallel configurations, 
in our numerical simulations we choose (V e , Vft)=(l, 1) and (V e , V^) = (— 1, —1), while for the 
antiparallel configurations we choose (V e , Vh)=(— 1, 1) and (V e , Vh)=(l, — 1). From our ex- 
tensive numerical results, we find that only one spin channel of EOPs is dominated for 
each of the four choices of (V e ,Vh). Furthermore, we find that the EOPs will obtain k- 
dependent phases due to the Rashba SO coupling. For convenience of discussion, we define 
Xfc=arctan(sin k y / sink x ) and Tfc=arctan(6fc/afc). As a typical example, the numerical results 
of EOPs for (V e ,Vh) = {l,l) are shown in Fig. 2. In this case, one can find from Fig. 2 
that the component A^(k) in EOP matrix Eq. ([6]) is dominant, while the amplitudes of 
the other three components (A^, A-j-jj., A^) are negligibly small. With keeping in mind that 
the fc-dependent phases of EOPs are obviously due to the Rashba SO interaction, we have 
analytically constructed various possible SO interaction-induced phases in EOPs and turned 
to compare these analytic approximate expressions with our exact numerical results. Table 
I summarizes the most optimal approximate phases for the four magnetic configurations. 
As an illustration, we plot in Fig. 3 our derived approximate condensate phases for the case 
of (V e , Vfc)=(l, 1), and compare them with the exact numerical result shown in Fig. 2. The 
agreement is clear. 

TABLE I. Dominant spin channel and approximate analytical EOP phase factors for different magnetization configurations 

(V e , Vh) dominant EOP component phases of (A-f^, A^, A^, A^) 
(1,1) A^(jfe) (l,-ie iTk ,ie iXk ,e i( - Xk+Tk ^) 

(-1,1) A n (k) (ie- iXk ,e- iXk+iTk ,-l,ie iTk ) 

(-1,-1) A n (k) (e-^ Xk+Tk \-ie- iXk ,ie- iTk ,l) 

(1,-1) A u {k) (-ie iTk ,-l,e iXk - iTk ,-ie iXk ) 



A n (k) A n (k) 
A u (k) A u (k) 



(6) 
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FIG. 2: (Color online) The left (middle, right) panels respectively show our calculated magnitudes 
(real parts, imaginary parts) of the EOPs A^(fc), A^(k), A^(k), and Ajjj.(/c), at a typical setup 
of magnetization parameters (y e , Vh)=(l, 1). 
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FIG. 3: (Color online) (from left to right) Phases factors of the EOPs A^(k), A^(k), and A^(k) 
that are listed in Table I with (V e ,Vh)=(l,l). The imaginary part of the component A^(k) is 
negligibly small (see Fig. 2(a2)). Here the upper and lower panels respectively plot the real and 
imaginary parts of these three phase factors. The black cirques denote Fermi surface. 
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With the help of Table I, we expect that the fc-dependent phases in the EOPs may 
lead to the nontrivially chiral topological orders. For instance, let us consider the case of 
(V e , V h ) = {l, 1). In the continuum limit, e iXfe ~^±^, and thus A lt (k)^i\A lt (k)\^^. That 
means the (p x +ip y )-like pairing emerges. 

Moreover, an explicit picture of chiral TEI can be well understood in the two-band 
approximation. To reveal this fact, first the non-interacting part in the total Hamiltonian 
is rewritten in the single-particle eigenstate space as 

n = J2 E es(k)4s(k)^es(k) + EUk^hsi-k^U-k), (7) 

k , s 

where E es =(f -/i e +s V « 2 (sin 2 k x + sin 2 k y )+V e 2 and E^-^+Hh+sy/ ' P\a\ + b 2 k )+V 2 
(s=+, — ) are respectively electron and heavy-hole band energies, and ip ps denotes the rele- 
vant annihilation field operators. Here the single-particle eigenstates are given by 

^(k) = & H/ + (fc)e-^,/_(A;),0,0] T , 
if e+ (k) = e"* [f4k),-tf + (k)e^,0,0] T , 
<p h -(k) = e id - [0,0,tg + (k)e^,g.(k)] T , 

<p h+ (k) = e**" [0A9-(k)J9+(k)e-^] T , (8) 



i , /n a-v/sin 2 fec+sin 2 k y , 

where J±{k) = — / ¥ — and 

^a 2 (sin 2 k x +sm 2 k y )+(y/a 2 (sm 2 fc^+sin 2 k y )+V 2 ±V e ) 2 



q+(k)= i v fc fc Note that 9u and du are fc-dependent phases 

^^al + bl) + ^P(ai+b 2 k)+ V 2 ±V h ) 2 

and are in principal determined, during exciton formation, by exactly solving the ground 
state of the system through our above self-consistent calculation. The single-particle bands 
E p ±(k) are shown in Fig. 1(b) (solid curves), from which it is easy to find that the excitons 
are preferably formed between the lower electron band E e _ and the upper hole band Eh + . 
Moreover, the pairing relates to the Fermi surface of the bilayer system. With the values 
of the tunable parameters shown in the caption of the Fig. 1, the band E e _ and band 
Eh + have the nearly perfect nesting Fermi surface with the Fermi energy Ep being nearly 
zero, (namely, fi p =—4t p ). Hence, we can deal with pairing in BCS picture in this situation. 
Now, we consider the electron-hole interaction part in Eq. (j2J) in terms of the filled electron 
band E e _ and hole band Eh+- In order to obtain an explicit picture, we use a rough 
approximation by assuming a short-range interaction potential U^ eh \q)=U5(q). Then, after 
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mean- field treatment, the resultant two-band Hamiltonian for our exciton system is given 
by 

H MF « ^S e _(fe)^_®Ve-® + ^E h+ (k)l() h+ (-k)lPl + (-k) 



- V(A(^)C(^)4 + (-^) + H.c.) + 



A(k) 



U 



(9) 



where A(k)=^^U f 2 (k)g 2 s ,(k) (^h+{—k)%l) e _(k)j (s,s'=±). The straightforward calculation 

s,s' 

can prove ^^/ S 2 (fc)^/(A;)~1 near the Fermi wave vector kp. So A(A;)^A is almost k- 

s,s' 

independent and only nonzero around kp in BCS-type picture. In practice we can introduce 



a factor r ){k)=e 



-c(k—kp) p iui 



e %UJ (c and oj are real constants) to fit our exact multi-band self- 



consistent numerical results (say, Fig. 2) in the whole BZ. The gaped energy spectrum of 
Hmf is shown in Fig. 1(b) (dashed lines). Now, in the two-band approximation, the EOPs 
in Eq. ([6]) have the expressions as follows: 



A(k) 



if+(k) ( k)e~ iXk f + (k)g + (k)e- iXk+iTk 
-f-(k)g-(k) if^{k)g + {k)e^ 



(10) 



where the phases e % ^ k are confirmed through our self-consistent calculation. It turns 

out that Equation (fTOil gives a nice description of the numerical results. 



III. CHIRAL TOPOLOGICAL ORDER 

In the presence of exchange fields (V e , V4) induced by the ferromagnetic films, the TRS of 
the system is broken. No less than the AQHE, the nonzero TKNN number can undoubtedly 
characterize the topological nature of the system if a stable bulk gap separates the ground 
state and excited states. That means the topological property of the system will not be 
changed without bulk gap closing in spite of adiabatically deforming |A CT(T /(fe)| at the given 
exchange fields. Hence, j(k) in Eq. ffTU]) is inessential for the system's topological property. 
Moreover, only the dominant component of EOPs decides the system's topological property 
at the given (V e , Vh,)- The straightforward calculation of Itknn in Eq. ( TTTT) can prove the 
above two arguments. 
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In the following discussion, we use Eq. (fTUj) to consider the system's topological prop- 
erties. In general, in the spin-dependent Nambu space (er f in 
different spin channels are affected by the effective exchange fields. Additionally, the strong 
Rashba SO interaction flaws the spin polarization of the carries along the z direction. The 
total effect leads the factors f±{k)g±(k) to emerge in different spin channels of EOPs, and 
which decide the dominant one at given (V e , Vh,)- For convenience of the following discussion, 
we use (A™, Af, A**, A d d ) to denote A (f + (k)g_(k), f + (k)g + (k), f_(k)g_(k), f_(k)g + (k)). 
The topological nature of the ground state \uo(k)) can be charactered by non-zero Itknn, 



which reads 

Itknn = ~ [ d 2 k{l (k), (11) 
2tt J bz 



where ft (A;)=-2Im < 



l^-y is the ground-state Berry curvature in BZ. The results are 
summarized in Table II, which definitely shows chiral topological order with its winding 
behavior depending on the choice of exchange-field parameters. From the bulk-edge corre- 
spondence, the nontrivial bulk topological number implies gapless edge states in the system 
with finite size. 

TABLE II. The TKNN numbers for effective exchange fields and corresponding EOP amplitudes 

(V e , Vh) A (Ag", Af, A du , A d d ) I TKNN 
(1,1) 0.5 0.5(0,0,1,0) 1 
(-1,1) 0.5 0.5(1,0,0,0) -1 
(-1,-1)0.5 0.5(0,1,0,0) -1 
(1,-1) 0.5 0.5(0,0,0,1) 1 

In order to confirm the existence of the gapless edge states, we assume that the square- 
lattice system has two edges in y direction and is boundless in x direction. Correspondingly, 
we choose open boundary condition in y direction and periodic boundary condition in x 
direction of the lattice Hamiltonian in Eq ([T]) in mean-field approximation. The calculated 
energy spectrum at a typical case of (V e , V/ l ) = (l, 1) is illustrated in Fig. 4 (a). The red-solid 
and blue-dashed lines correspond to the different edge states with contrary chirality. It is 
easy to find that the number of the gapless edge states is consistent with the bulk theory 
characterized by Itknn- 
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FIG. 4: (Color online) (a) The energy spectrum of the bilayer square-lattice system with two edges 
at the y direction. k x denotes the momentum in the x direction. The magnetization parameters 
are set at (V e , Vh)=(l, 1). The red-solid and blue-dashed lines denote the edge states locating on 
different edges, (b), Six-terminal Hall bar for detection of the edge states. The red-solid and 
blue-dashed lines with arrows represent the edge modes propagating in opposite direction. 

IV. TRANSPORT PROPERTY OF EDGE STATES 

The nontrivial transport phenomena can be predicted due to the emergence of the edge 
states in our system. From Fig. 4 (a), we can find that the edge sates in different chiral 
topological order propagate on each boundary with opposite velocities and can be described 
by the following Hamiltonian 



where ± represents different edges and 77=1,..., 4 labels four different kinds of magnetic config- 
urations, namely, Ai=A4=l and A2=A3=— 1, vp is the Fermi velocity and k x is the momentum 
measured from the Fermi surface. The quasiparticle operators for case (V e , V/ l )=(l, 1) read 




(12) 



k x >0 



Ji{k x , y) = u\{k X i y)et(y) + v x {k x , y)h\{y). 



(13) 
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The other cases have the similar forms. Due to the missing of the PHS, the quasiparticles 
are not Major ana fermions. 

The edge states in the AQHE systems can be usually detected through the Hall conduc- 



tance responding to the external electromagnetic field [16| 17|. However, the edge states 
in our system are excitons which are charge neutral. A simple approach is to use thermal 
transport measurement which is often used to judge the pair properties in high-T c supercon- 



ductors 



181 ] 191 ] . The six-terminal Hall bar showed in Fig. 4 (b) for detecting the edge states 



of quantum (spin) Hall effect can be used to detect the thermal conductance. The same 



setup has been used by Sato et. al 



201 ] to detect the edge state in topological superconductor. 



We give the similar considerations with that in Ref. 20| as follows. The temperature must 
be sufficiently lower than the exciton gap (T-cA ) in order to suppress the contributions 
from the fermionic excitations (electrons and holes) in the bulk and bosonic (phonons) ex- 
citations. The thermal conductance is defined by G(T)=J 14 (T)/(AT) 14 , where Uj{T) is a 
thermal current between contacts % and j, and (AT)ij is the temperature difference between 
these contacts. In the low temperature limit, the T-dependence of G(T) have three origins: 
the linear law ocT from edge states for phase r/, the exponentially low ~e~ A °/ T from bulk 
quasiparticles and the power law ocT 3 from phonons. Furthermore, in analogy with the 
quantum spin Hall current discussed in Ref. 21[ , there is no temperature difference between 
contacts 2 and 3 (5 and 6) because the edge current is dissipationless. 



V. CONCLUSION 



In conclusion, we have presented a scheme to realize the chiral topological excitonic insu- 
lator in the double quantum wells adjacent to two ferromagnetic films. We have predicted 
different topologically nontrivial orders emergent along with changes in the magnetization 
orientations in the ferromagnetic films. The topologically nontrivial orders can be charac- 
terized by the chiral topological numbers defined with TKNN numbers in bulk system or 
chiral edge states in edged system. Furthermore, we have given an experimental scheme to 
detect the excitonic gapless edge states. 
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